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GEOMAGNETIC RADIOPLASi! UNFOLL <GRUF) 


by 


John S. Malik 


ABSTRACT 


A method of invjiting the geomagnetic component of the 
radioflash oignai from a nuclear explosion to obtain the 
gamma-ray timo history was proposed by K. 0. Dracott of 
the Atomic Weapons Research Establishment. A simplified 
development of an elaboration by 8. R. Suydan has been 
programmed for small calculators in a Corn suitable for 
interim field analysis of such data. The development of 
the program is contained in the report. 


1. 
INTRODUCTION 


The most useful component of the radio- 


flash signal emitted during a nuclear weapon 


explosion for diagnostic purposes in the geo- 


magnetic component. As its name implies it 


results from the deflection of the Ccmpton 


electrons, produced by interaction of gamma 


rays with air molecules, by the static geo- 


magnetic field of the earth. It has the 


character of a magnetic dipole signal; the 


radiated electric field is perpendicular to 


the geomagnetic field. Its generation has 


been described in several reports,notably by 


Suydam. 
Its potential for the determina- 


tion i-J the gamma time history of the source 


producing the signal was recognized first by 


E. D. Oracott.2 The implementation of the 


proposal has been by Moody and Hill with 


further extension by Suydam. 
Suydam's meth- 


od has been implemented by Malik in the Los 


Alamos Scientific Laboratory (LASL) program 


HMSD. Program GRUF is a simplified unfold- 


ing technique which has been programmed for 


the Wang 700 and Hewlett-Packard 9820 calcu- 


lators; it is suitable for early data reduc- 


tion while in the field. As in the develop- 


ment of the program HMSD, the development 


of the program GRUF closely follows 


Suydam's. 


II. DEVELOPMENT 


The starting point is Maxwell's equa- 
tions in MKS units: 
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In polar coordinates, the magnetic or TE set 


of equations is 
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Transforming to retarded time coordinates: 
r ' = r, T = t-r/c, for which the operators 
become j^ •• ^ , - — ^ >•§•£•* yr and dropping 
the prime on r, the Bo equation becomes: 
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The B Q and 9 r £qa. f«) and (Si are inte- 
grable to 
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In the transformed coordinates, C<j. (3) be- 
comas 
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Substituting for 9 r and! B@ in Eq. <9) 
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At larg* distances {o -OS ,:•;•• 0 and 
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Recognizing th* asimuthal dependence, we 
will omit it temporarily and solve for 
r£. it i • 17/2. The gamma-ray dependent 
funcelone are approximated by 


(17) 


By restricting the solutions to shovt times, 
the terns involving integrals mr.y be t'ropped: 
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This is the starting point for Suydam's de- 
velopment . 
Choosing the polar axis antiparallel to 
the earth's magnetic field, the transverse 
current J. at the azimuthal angle 0 is given 
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where u is the electron (nobility, W^ is 
energy per ion pair, e the electron charge, 
and *e, \m are y-ray Compton ard absorption 
mean free paths. The function g has dim«n~ 
•ion* of energy; f, that of energy/tine. 
They are connected by the air-ion equation: 


where t is the range of the Compton electron, 
a the Larmor radius, Jc the Compton current. 
Equation (11) integrates to 
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whertt a is the electron attachment rate to 


0,: n it; th<- only important loss rate for 


tin- t ttwa "f interest. The electron mobil- 


ity i;(r »t; it la ted to conductivity by « c neu > 
with n beinq the electron density. 


Denoting time derivativea by dot, E. as 


defined by Siq. (141 btjccir.os 


where v is the number of secondary electrons 


produced per primary v-ray of energy W. . 


Putting Eqs. (26) and (27) into Eq. (21), 
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kritinq the equation for rt:,, Kq. (id as 
I s- • >: (In ! I (1+Qn 
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where S,. it; the value of r satisfying 
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The correction quantity 0 
is tabulated by 
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Equation (23) in the form 
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gives the differential equation for x or R s. 


Using Eq. [7.1) gives 9 ( 1 ) . Combining Eqs. 
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and '„ t.1; the y-ray transport mean free path. 


By logarithmic differentiation: 
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The ratio V/:. is given by 
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Prom a measurement of xZ , integration 
of Eq. (32) followed by a substitution into 
Sq. (34) yields tho gamma-ray tine depen- 
dence f. This assumes a knowledge of the 
other parameters in these relations, most of 
which are functions of the Y-ray energy. 
Since the source covers a spectrum of Y-ray 
energies from below 0.1 HeV to perhaps 
10 MeV, a weighted average over the spectrum 
is needed. As the variable x - R/*T is not 
a linear function of distance, the averaged 
variables were fitted as functions of x 
rather than the radial distance. The elec- 
tron mobility is approximated by 


Eqs. (32) and (34) by 


rE^ » REV/ (cos* . sine) 


The angle 0 is the supplement of the true 
bearing to the source minus the magnetic 
declination. 
The integration is performed by a sec- 
ond order Runge-Kutta scheme: 
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where p is the molecular per cent of water 
in the air. The electron attachment rate in 
approximated by 


S - 0.09 (•£-\na~1 


The correction factor for the integral is 
fit by 


Q^ « <S.O7 + 4.61x)/(l*16x+8.36x2) . 


The quantify K, as the result of the averag- 


ing is 
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The measured vertical electric field, E y, 
obtained at a distance R in meters co-alti- 
tude with the source at an aximuthal angle 8 
a:id at a region where the dip angle is i>, re- 
lates to the electric field rE. in the 
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or in this instance: 
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